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Abstract 

We report experiments on defect-tracking in the state of undulation chaos observed in thermal 
convection of an inclined fluid layer. We characterize the ensemble of defect trajectories according 
to their velocities, relative positions, diffusion, and gain and loss rates. In particular, the defects 
exhibit incidents of rapid transverse motion which result in power law distributions for a number 
of quantitative measures. We examine connections between this behavior and Levy flights and 
anomalous diffusion. In addition, we describe time-reversal and system size invariance for defect 
creation and annihilation rates. 
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Topological defects within patterns are observed in many systems to move 
in a spatiotemporally chaotic fashion. We examine the motion of such defects 
within a defect-turbulent state observed in thermal convection of an inclined 
fluid layer. We characterize the trajectories of the defects both by analogy to 
fluid turbulence — velocity distributions, diffusion, and power spectra — and by 
properties dependent upon the topological characteristics of the defects: pair 
creation/annihilation and interactions. 



I. INTRODUCTION 

Nonequilibrium systems with similar symmetries often form patterns which appear to be 
universal in spite of having been formed by different physical mechanisms |I| . In particular, 
reduced descriptions of the patterns often quantify the similarities in behavior so that under- 
standing of one system can lead to insights in multiple systems. A class of spatiotemporally 
chaotic states exhibiting defect-mediated turbulence || has been found in such diverse sys- 
tems as wind-driven sand, electroconvection in liquid crystals |§, nonlinear optics ||, fluid 
convection || ||, and autocatalytic chemical reactions J7|. In many cases, such systems have 
been modeled via the complex Ginzburg-Landau equation || ^ |10|, [TT], These various 



defect turbulent patterns are characterized by an underlying striped state which contains 
dislocations (point defects) where the stripes dead-end within the pattern. Locally, the de- 
fects distort the orientation and wavenumber of the stripes and the nucleation, motion, and 
annihilation of the defects constitute a spatiotemporally chaotic system. An example from 
inclined layer convection is shown in Fig. [I]. 

Previous work on defect turbulence has focused both on snapshots of such patterns |3|, £|, 



TU| and the dynamics and interaction || TT|| . However, there are numerous open questions 



about defect turbulence: characterization of the defect motions, interactions between the 
defects, and the extent to which the analogy to turbulence is appropriate. The degree to 
which such characteristics are similar in different defect-turbulent systems remains to be 
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FIG. 1: Sample Fourier-filtered shadowgraph image of inclined layer convection at e = 0.08 and 
7 = 30°. Black circle encloses a positive defect; white, a negative. Arrow indicates tearing region 
of low-amplitude convection. Uphill direction is at left side of page. Region shown is the subregion 
of size 5 lei x 63d used for analysis. 
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FIG. 2: Schematic of inclined layer convection with associated coordinate system. AT = T^ot — 

Tcold- 

explored. 

Investigations of pattern formation in variants of Rayleigh-Benard convection (RBC) have 
been particularly fruitful [[b| . The state of undulation chaos (shown in Fig. |l] and described 



in 



14| , |i~5| , |T6|1 ) observed in inclined layer convection (Fig. ||) exhibits defect turbulence and 
is well suited to investigations on the dynamics of defects since spatially extended systems 
and fast time scales are experimentally accessible. This allows for tracking of point defects 
through their creation, motion, and annihilation. In the observed pattern, the stripes contain 
undulations as well as defects; both are spatiotemporally chaotic (further characterization 
to be published in [|i~6f). 

A number of features stand out in the defect trajectories we observe in undulation chaos. 
As exemplified in Fig. |3|, the inherent anisotropy (due to the inclination) is apparent, with 
the trajectory meandering preferentially in the transverse direction. Occasionally, there 
is a rapid burst of transverse motion as marked by the black diamonds, corresponding to 
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FIG. 3: (a) Example of a long trajectory for a defect with positive topological charge at e = 0.08. 
Positive y indicates a downslope direction (see Fig. |l]). (b) Corresponding transverse velocity. The 
symbols ♦ mark the ends of a flight with displacement AX and duration At. Dotted line is the 
trigger velocity for detecting flights. 

a tearing of the pattern across the rolls. Such behavior appears be related to Levy flights 
fP7| ) |T^| , PHH for which distributions of step sizes display power laws. Furthermore, we are able 
to study the ensemble of trajectories to gain insight into defect pair creation, interaction, 
and annihilation. 



II. INCLINED LAYER CONVECTION 

In inclined layer convection (ILC), a thin fluid layer heated from one side and cooled from 
the other is tilted by an angle 7; the system is anisotropic due to the resulting shear flow (see 
Fig. |2|). The fluid becomes unstable above a critical temperature difference AT C . At fixed 
7, we describe how far the system is above the onset of convection via the non-dimensional 
driving parameter e = A ^ — 1. 

At low angles of inclination, buoyancy provides the primary instability (analogous to 
RBC) and the convection rolls are aligned with the shear flow (longitudinal) both at onset 
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and for the secondary transition to undulation chaos [p0|| . Over a range of intermediate 
angles (15° < 7 < 70° for Prandtl number a ~ 1), the longitudinal rolls become unstable to 
undulation chaos above e ~ 0.01. It is this defect turbulent state which we investigate; some 
of its properties have been previously described in other work [HI ffB]. Above e ~ 0.1, there 
is a further transition to a state of competing ordered undulations and undulation chaos. 
We examine trajectories from both sides of this transition. 

The apparatus used in this experiment is of the type described in de Bruyn et al. pi] , 
modified to allow for inclination. The fluid used was CO2 at a pressure of (56.5 ± 0.01) bar 
regulated to ±0.005 bar with a mean temperature of (28 ± 0.05)°C regulated to ±0.0003°C. 
As determined from a materials properties program pl| , the Prandtl number was o = 
1.140 ± 0.001. A cell of height d = (388 ± 2)/im and dimensions 203d x lOOrf was used, for 
which the vertical diffusion time was t v = d 2 /n = (1.532 ± 0.015) sec. The fluid was weakly 
non-Boussinesq conditions: Q = 0.8, as described in Jl3| for horizontal fluid layers. All 
experiments were performed at a fixed inclination of 7 = 30°, within the regime of buoyancy- 
instability. Images of the convection pattern were collected using a digital CCD camera, via 
the usual shadowgraph technique [f2T. 22]. Images were collected at 3 frames per second in 
one of two formats. Six- hour (2 x 10 4 t„, 80000 frames) continuous runs of data were obtained 
at two values of e: 0.08 (four runs) and 0.17 (two runs). For 17 values of e between 0.04 and 
0.22, short runs with 100 images were collected, separated by at least IOO77 for statistical 
independence. At each of these e, at least 400 repeats (up to 600 for the lowest values of 
e) were performed. Each value of e was reached by a quasistatic temperature increase from 
below. In addition, a run with quasistatic temperature decreases was performed between 
e = 0.12 and e = 0.06 to check for hysteresis, which was not observed. Only data from the 
central homogeneous region of dimension 51c? x 63<i was utilized during the analysis unless 



noted otherwise; see |16| for details on the choice of this region. Size-dependent effects of 
this region are discussed here. 



III. DEFECT TRAJECTORIES 



Here we consider only the defects themselves — an ensemble of moving, interacting, 
charged "particles" - - and not the underlying pattern. By reducing the data to a set of 
defect trajectories of known charge, birth, and death we can observe how these particles move 
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individually and as an ensemble. Single defects may enter or leave through the edges of the 
system, or be created/annihilated in pairwise events between defects of opposite topological 
charge [TJ|. During their lifetimes, the defects have positions and velocities which may be 



related both to the underlying pattern and to the presence of other defects. We examine 
properties averaged over these effects, and also isolate some effects due to the latter. 

The topological defects are located at points where a roll pair ends within the pattern 
(see Fig. [I]). At these points, there is a phase discontinuity and the topological charge is 
determined from the phase jump along a contour around the defect: 

V0 • ds = n2n, (1) 

The phase field <f>(x, y) is determined from Fourier-demodulated images for which the com- 
plex field ip(x,y) = \A(x,y)\ e % ^ x,y > has been reconstructed. We detect the locations of 
defects by finding points where 9ft(?/>) = and = and a phase discontinuity occurs. 

(Details of this technique are described in |TB|.) In undulation chaos, we observe only defects 
with n = ±1; examples of each topological charge are marked by circles in Fig. [I]. 

After the defects were detected in each image, the defect locations were connected to 
form trajectories. This was done by matching each defect already on a trajectory to its 
closest same-signed neighbor in the subsequent frame. If the reverse process also agreed 
on the match, then the subsequent defect was added to the list of defects located on that 
trajectory. Similarly, the ends of positive and negative defect trajectories were matched 
up to obtain the locations of all pair creations and annihilations. When the analysis was 
completed, broken trajectories (those missing a creation or annihilation) were rejoined. 

Since the defect positions were determined only to the nearest pixel, the data was later 
smoothed to interpolate the coordinates along the trajectory. Each X{ = x(tj) and yi = 
y(ti) along a trajectory was replaced with the weighted average of its neighbors, using the 
Gaussian weighting function 

W *i = W (2,)3/ 2 eX P[ 25T"] 

■exp[-^] -exp[-^] (2) 

This resulted in a smoothed trajectory with coordinates Xi = w ij x j an d fji = ^2 w ijVj 
(tildes will be dropped hereafter for convenience). The appropriate widths o in Wij were 
determined by examining a range of parameters and finding convergence for a x = a v = 2 d 




FIG. 4: Velocity autocorrelation in (a) transverse and (b) longitudinal directions at e = 0.08 and 
e = 0.17. Data is shown for positive defects; results for negative defect trajectories were similar. 

and at = 2r v . The advantage of this weighting method is that it automatically adjusts 
the fit length so that along faster-moving segments of the trajectory we average fewer data 
points than along slower-moving segments. 

From the trajectories x(t) and y(t), we applied a local weighted linear fit (again using Wij) 
at each point along the trajectory to obtain the corresponding velocities v x (t) and v y (t). We 
also determined v = v 2 + v 2 and 6 = arctan ^, corresponding to a downslope (y) direction 
at 9 = 0°. A raw trajectory is shown in Fig. |3], along with the smoothed trajectory and the 
transverse velocity component v x (t). 

Defect turbulence in inclined layer convection is anisotropic, with v x corresponding to 
glide motion across the rolls and v y to climb motion along the rolls. These are fundamen- 
tally different behaviors: the former adjusts the orientation of the rolls, and the later the 
wavenumber. As expected, the decomposed motions are found to be poorly correlated with 
each other, with a linear correlation coefficient R 2 = 0.16 for positive defect trajectories and 
R 2 = 0.12 for negative. Therefore, we separate all defect positions and motions into their 
two components and examine each direction independently. 

Defect velocities are observed to be correlated over time scales of less than « 10r„, as 
determined from temporal autocorrelation functions of v x (t) and v y (t). The autocorrelations 
are plotted in Fig. [|, averaged over the ensemble of trajectories. The defect motions exhibit 
short-range order, with exponential decay during the time before the zero-crossing. 
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IV. DEFECT MOTION 



We characterize the meanderings of the defect via its central moments, represented by 



integer powers of the deviations from its mean value. For trajectories Xi(t) with mean Xi(t), 
the n th moment is given by the quantity 

^t) = ((x t (t)-xi) n ) (3) 

where (•) represents an ensemble average, and a similar calculation can be made for the y(t). 
For normal diffusion or random walks, fj, n oc t n ^ 2 , where [12 — 2-Dt provides a value D for 
the diffusion constant. 



Investigations in a broad variety of situations [see, for example, |17], [Tg, [L9[ and 
references therein] have found anomalous diffusion instead, where LI2 = 2Dt a . For a < 1 a 
system is said to be subdiffusive, a > 1 is superdiffusive, and a = 2 is ballistic transport. 
Anomalous diffusion is frequently associated with Levy walks/flights, in which the PDFs of 
the lengths and durations of the flights are power laws and have infinite moments. Since the 
central limit theorem no longer holds in such cases, the probability of long jumps (waits) 
will enhance (retard) the diffusive behavior. Similarly, the power spectrum S(u) of the 
associated velocities will exhibit power law behavior as well. Some recent examples of such 



behavior in fluid dynamics include tokamak density fluctuations [2^ and particle motion in 
two-dimensional rotating fluid flows |25|. Defect random walks have also been considered 
for models of intracellular Ca 2+ dynamics ||, in which subdiffusion was observed. 

In Fig. |5| and |^ we calculate moments for those values of t at which 100 or more trajectories 
have at least that duration. We look for behavior of the type fi2 oc t a , 114 oc t 13 , and liq oc t 7 
in both the x and y directions. In the variance (^2), superdiffusive behavior is evident at 
both e, particularly in x where a ~ 1.5. The exponents for y are lower but also indicate 
super diffusion. The exponents of 114 and fie do not suggest super diffusion in x, and the y 
behavior does not follow a power law. 

We locate defects flights as the portions of the trajectories which occur while \v x \ > 1 d/r v . 
An example flight is shown in Fig. || For each of these flights, we determine the duration 
At and displacement AX. Distributions of these values are shown in Fig. [7| Due to the 
constraints of a finite system size, the maximum observable AX is limited. Nonetheless, 
the distribution of flight times At is observed to have behavior consistent with a power 
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FIG. 5: Displacement moments (Eqn. ||) and local fits to the associated exponents in x (black) 
and y (gray) directions at e = 0.08. Dashed lines are values for normal diffusion or random walks, 
[i n (x t n l 2 . Data is shown for positive defects; results for negative defects were similar. 

law tail of exponent approximately —3. For AX, behavior suggestive of a power law with 
exponent approximately —4 is observed, indicating that there is not a characteristic velocity 
associated with the flights. 

Physically, these flight events correspond to a tearing of the rolls when the pattern shifts 
along a line in the x direction, so that each roll broken in the process rejoins with the one 
next to it. Such a tearing line is visible next to the white arrow in Fig. [|, where there is 
weaker convection along the slopes of the undulations. Since each roll must move by only 
2d to reconnect with the next roll, defects can be transported rapidly from one end of the 
tear line to the other. 

By analogy to fluid turbulence, we plot probability distribution functions (PDFs) and 
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FIG. 6: Displacement moments (Eqn. ||) and local fits to the associated exponents in x (black) 
and y (gray) directions at e = 0.17. Dashed lines are values for normal diffusion or random walks, 
[i n (x t n l 2 . Data is shown for positive defects; results for negative defects were similar. 

power spectra of the defect velocities. Figs. |8| and |] show the results for distributions in x 
and y. The PDFs of v x are independent of topological charge and show steep power law 
tails with an exponent of approximately —3.5. This is consistent with the results obtained 
for defect flights displacements and durations. As can be seen in Figs. || and |], the peak 
of the distributions have a Gaussian shape, indicative of smoothing due to noise. Each 
distribution is also of a shape described by Tsallis statistics of the type described in |^6[ . The 
PDFs for v y show lower characteristic velocities and (particularly at e = 0.08) dependence 
on the topological charge. Fig. ID shows a complementary picture in (v, 9) space via a 
two dimensional histogram. The x flights are now visible as the ridges at approximately 
9 = ±90°. This transverse direction not only comprises most of the defect motion, but is 
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FIG. 7: Distributions of flight (a,c) durations At and (b,d) displacements AX at e = 0.08 and e = 
0.17. Positive defects are +, negative defects are squares, black points are from entire trajectories, 
and gray from segments within the homogeneous subregion. Dotted lines represent fit region for 
power laws. 



also the direction along which motion is the fastest. 

A second effect visible in Figs. [|, ^, and [10] is an asymmetry in the ±y behavior which 
shifts the (v,6) peaks in the upslope or downslope direction, away from ±90°. Although 
the e = 0.17 defects are more isotropic in their motion than the e = 0.08 defects, in both 
cases the negative defects prefer the downslope direction over the upslope, and the positive 
defects prefer the upslope. In fact, Figs. |] and ^| show that this asymmetry is the dominant 
behavior at lower e: it is rare for defects to travel other than in the preferred direction. 
Asymmetries in defect behavior with respect to topological charge have been observed in 
other anisotropic systems such as Langmuir circulations f27|, although for quantities such 
as the surface convergence rather than the defect velocity. In this case, the asymmetry may 
be related to the breaking of the z symmetry by non-Boussinesq effects. 

In a system consisting of Levy flight behavior, the power spectra S(uj) of v(t) is expected 
to show power law behavior for 7r 3> uo ^> uo* = 1/At*, where At* is a short time associated 



with the low-t limit of power law behavior of the PDF of At for flights [19]. Such a cutoff is 
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FIG. 8: (a,b) Longitudinal and (c,d) transverse velocity PDFs at e = 0.08 on log-log and semi-log 
(for v > 0) axes. Positive defects are shown with black + and negative defects with gray squares. 
Approximately 10 6 defect velocity measurements were used to determine the probabilities in each 
graph. Dotted lines are Gaussian fits to whole data, with standard deviations a y = 0.08 and 
a x = 0.59. Dashed lines represent fit region for power laws. 



relatively large (~ 5r„) in the defect flight durations shown in Fig. 0, limiting the range over 



which we would expect to see power laws in 5" (a;) for defect trajectories. Fig. [TT] suggests 
random- walk behavior (constant, uniform spectrum) for low u, crossing over to a flight- 
related power law at high u>. The exponent in all cases (e, x, y, and topological charges) is 
approximately —1.8. 
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FIG. 9: (a,b) Longitudinal and (c,d) transverse velocity PDFs at e = 0.17 on log-log and semi-log 
(for v > 0) axes. Positive defects are shown with black + and negative defects with gray squares. 
Approximately 10 6 defect velocity measurements were used to determine the probabilities in each 
graph. Dotted lines are Gaussian fits to whole data, with standard deviations a y = 0.10 and 
a x = 0.51. Dashed lines represent fit region for power laws. 
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FIG. 10: Histogram for (v, 9) pairs. Negative defects plotted with negative speed and positive 
defects plotted with positive speed. Grayscale is shown on a logarithmic scale, (a) e = 0.08 (b) 
e = 0.17. 
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FIG. 11: Power spectra of (a, c) v x and (b, d) v y at e = 0.08 and 0.17 for positive (black solid) 
and negative (gray dashed) defects. Dotted lines show power law fit region for positive defects. 
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FIG. 12: Schematic diagram of determination of defect pair separations in y, with associated Ay 
measurements. 
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FIG. 13: PDF of defect pair separations (Ax and Ay) for strips of width 5y = 5x = 2d as shown 
in Fig. [T^. (a,c) Transverse separations and (b,d) longitudinal separations. Solid lines are (+, — ) 
pairs, dotted are (+,+), and dashed are (—,—). Approximately 10 5 elements used for each curve. 

V. DEFECT INTERACTIONS 



The defects carry topological charge, and the resulting phase fields allow them to interact 
with each other. By examining the relative position of nearby defects, we can gain insight 
into their interactions. We examine strips of width Sy = 2d (or 5x) and calculate the 



distances Ax (or Ay) separating all pairs of defects within the strip (see Fig. [12| for a 
schematic). 



15 



Fig. [13] plots PDFs of this data for all (+,—), (+,+), and (— , — ) pairs of defects: 
V(Ax | 5y < 2d) and V(Ay\5x < 2d). The data is limited by the finite size of the re- 
gion under consideration, and is only plotted over half the dimensions of the region. For 
both e, it is rare for oppositely-charged defects to be located within 3d of each other. This 
is associated with the annihilation of such pairs which pass too close in either direction. 

Defect pairs of (+, +) or (— , — ) charge behave in similar fashions to each other, demon- 



strated by the closely overlapping dashed and dotted lines of Fig. [I3|. There is a repulsive 
effect in the x direction, with Ax increasing away from 0. The turnover in the e = 0.17 
PDF may be due to finite size considerations. In the y direction, there is again an excluded 
region of around 5d, an effect which is strong enough that secondary and tertiary peaks at 
lOd and 15d are visible as well for e = 0.17. 



VI. CREATION AND ANNIHILATION EVENTS 

A general feature of defect turbulent systems is the production and loss of defects through 
two mechanisms: pair creations/annihilations and flux through the boundaries. A mean- 



field approximation has been postulated by Gil et al. |T0[ and experimentally examined 



by Falcke et al. || and Daniels and Bodenschatz |15[ in which these rates depend on the 
number of defects in the system. 



C(N) 


= C 


(4) 


E(N) 


= E 


(5) 


L(N) 


= L N 


(6) 


A(N) 


= A N 2 


(7) 



where C(N) is the probability of a pair creation event happening per unit time, A(N) is 
pair annihilation, L(N) is a single defect leaving, and E(N) is a single defect entering. N 
is taken to be either the number of positive defects iV + or the number of negative defects 
N_, quantities which are approximately equal on average. Fig. [TJ] shows PDFs of the net 
charge (iV + — AL) for the observed images, with a mean of zero (topologically neutral), 
independent of e. The constant creation and entering rates can be understood as being 
generated by random events, independent of the number of defects already in the system. 
The annihilation rate scales as the number of positive defects times the number of negative 
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FIG. 14: PDF of net charge (N+ iV_) found in subregion for nine values of e. Solid line is a 

Gaussian fit. 

defects. The leaving rate is proportional to the number of defects present in the system. 

Based on this approach, we earlier derived and tested a universal distribution for iV which 
agrees with the experimental findings [|15|. We find, as well, that the subscripted coefficients 
depend on e, system size, and other physical parameters in systematic ways. 

The stationary distribution for iV was found using a recursive relation 

\oss(N)V(N) = gain(iV- 1)V(N - 1) (8) 



to describe the probabilities at adjacent iV [pL5| . In fact, the more stringent condition of 



detailed balance holds as well. By time-reversing all trajectories, creations become annihi- 
lations and entering defects leave the subregion. Analysis of the reversed trajectories for 



e = 0.07 is shown in Fig. |T| in comparison with the original data and found to exhibit the 
same behavior. 

Such defect gain and loss rates properly scale with the size of the region under consid- 
eration. We define a set of size-independent coefficients c , a , e , and l to quantify this 
behavior. Both creation and entering are random events, which have some rate per unit 
area and length, respectively. Thus, cq = and e = where S and P are the surface 
area and perimeter of the subregion, respectively. With primed and unprimed variables 
representing different subregion sizes, c = = and e = = For two subregions 
with the same density of defects n = N/ S, the annihilations per unit area (at constant n) 
is also constant. For A(N) = A N 2 = A n 2 S 2 = A' n 2 S' 2 and constant density of defects n, 
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FIG. 15: Comparison of gain/loss coefficients for forward-time (triangles) and backward-time 
(squares) trajectories at e = 0.07. Lines are fits to the expected ^-dependence as described 
in text. 









! c — 


— v 

— 


r 






— □ . 


■□— 


— B E 






A— 








G— 


— e € 


i A — 








i e — 


— 



15 20 25 30 
S (10 2 d 2 ) 




160 180 200 220 

P(d) 



15 

ra°10 
5 





15 20 25 30 
S (10 2 d 2 ) 



2.5 
2.0 
,1.5 
1.0 
0.5 
0.0 



o — e — $ — e — o 



a — a- 



160 180 200 220 
P(d) 



FIG. 16: Size-independent gain/loss coefficients vs. size of test region, e from bottom to top are: 
0.06 (O), 0.09 (A), 0.12 (□), 0.15 (0), 0.18 (y), 0.22 (x). Data is shown for positive defects; 
results for negative defects were similar. 



ao = oc A S. Analogously, l = L Q P. 

Fig. |l^ compares these rescaled coefficients for various sizes of test regions. Because of 
the x-y anisotropy, regions with the same aspect ratio as the homogeneous subregion were 
used. All four coefficients show constant rescaled rates for a given e. There is, however, 
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FIG. 17: Rate coefficients as a function of e for homogeneous subregion. A represent data taken 
under quasistatic temperature increase; V represent data taken under quasistatic temperature 
decrease. 

a slight trend in the annihilation rates: larger boxes allow more defects to find each other 
and annihilate, particularly at higher e. This may also relate to the larger test regions 
approaching the less homogeneous sidewall regions (see fll6l ). 

Fig. [T7| shows the coefficients Co, Eq, Lo, and Aq as a function of e. The error bars were 
determined via a bootstrap method: by re-sampling the data with replacement we obtained 
a distribution of values for each coefficient from which to estimate the error. At present 
we have no explanation for the particular shapes of these graphs, which appear to be linear 
to first order. Both Co(e) and Eo(e) appear to intercept the e-axis close to the e c ~ 0.02 
onset of undulation chaos |y|. For A (e) and L (e), both lines extrapolate to the origin, 
the expected behavior for N = 0. The data at e = 0.04 and 0.05 have been disregarded due 
to finite-size effects which are significant at low values of e since the undulations are a long- 
wavelength (k = 0) instability. Defects are increasingly rare for e — > e c as the wavelength of 
the undulations becomes on the order of the size of the convection cell. 



VII. CONCLUSION 

The defect trajectories observed in inclined layer convection display many intriguing 
behaviors. The occasional rapid motion of defects across many convection rolls can be 
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associated with power laws in various quantitative measures of their motility: diffusion, 
velocity PDFs, flight size PDFs, and velocity power spectra. While the system is strongly 
anisotropic and its defect motions are dominated by transverse motion, PDFs of the defect 
separations reveal significant correlations in the longitudinal direction as well. Finally, the 
trajectories allow us to probe the gain and loss of defects in which much of this interesting 
behavior can be averaged over to reveal general results. 

We have focused primarily on characterizing the observed flight behavior. If these rapid 
motions are, in fact, related to Levy flights then the exponents found in each of the power 
laws should be related to each other as has been presented in the literature for various 
random walk formalisms [18, 19, 28fl . However, comparisons to existing theory are difficult 
due to both the strongly anisotropic nature and the finite size of this system. Further 
investigation into this phenomenon, possibly in other systems, will undoubtedly be fruitful. 

The trajectories analyzed here were obtained for e = 0.08 and 0.17, of which the former 
is in a state of undulation chaos and the later intermittently exhibit ordered undulations. 
The dynamics of this transition is discussed in fl6| , and the defect trajectories should also 
be examined in light of the existence of the two states. The statistics presented here will 
presumably vary according to whether the system is in an ordered or disordered state. 
Furthermore, the motion of the defects during the transition may shed light on the nature 
of the transition. 

Other questions remain regarding the relation between the defect motion and the under- 
lying undulation pattern. For instance, the relationship between the motion of the defects 
and the stability of the local wavenumber. A related issue is that defects have been observed 
to "bounce" off regions of ordered undulations, sharply reversing direction. 

Finally, the relationship between the relative position and relative velocity of defects has 
not yet been investigated. Further work in this area would provide information about the 
attraction and repulsion of defects, particularly near creation and annihilation events. 

We are grateful to G. M. Zaslavksy, T. Geisel, J. P. Sethna, C. L. Henley, C. Beck, and 
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